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ABSTRACT 

A distinguishing trait of the three known Galactic recurrent novae with the shortest orbital 
periods, T Pyx, IM Nor, and Cl Aql, is that their optical decline time-scales are signihcantly 
longer than those of the other recurrent systems. On the other hand, some estimates of the mass 
of the ejecta, the velocity of the ejecta, and the duration of the soft X-rays emission of these 
systems are of the order of those of the other recurrent systems and the fast classical novae. 
We put forth a tentative explanation of this phenomenon. We propose that in these systems part 
of the material transferred from the companion during the hrst few days of the eruption remains 
within the Roche lobe of the white dwarf, preventing the radiation from ionizing the ejecta of 
the system and increasing the optical decline time-scale. We explain why this phenomenon 
is more likely in systems with a high mass transfer rate and a short orbital period. Finally, 
we present a schematic model that shows that the material transferred from the companion is 
sufficient to absorb the radiation from the white dwarf in these systems, ultimately supporting 
this scenario as quantitatively realistic. 

Key words: accretion, accretion discs - novae, cataclysmic variables - stars: individual (T 
Pyx, IM Nor, Cl Aql) 


1 INTRODUCTION 


2 T PYX, IM NOR, Cl AQL, AND THE OTHER NOVAE 


Classical novae are the result of a thermonuclear runaway induced 
in the envelope of a mass-accreting white dwarf. The white dwarf 
is not destroyed by the process and the system may undergo subse¬ 
quent eruptions. Recurrent novae, the rarest subclass, are those sys¬ 
tems that have been observed to experience more than one historical 
eruption, with a recurrence time-scale of a few decades. These are 
characterised by a high mass transfer rate as well as a high white 
dwarf mass (see e.g. |Anupama|2008^ . Only ten Galactic recurrent 
novae have been directly observed during at least two of their erup¬ 
tions; the statistics of these observations are reported by [Schaefer] 
( |2010| > and are partially summarised here in table 

In this paper we highlight an interesting feature of the data in 
table[T] the three recurrent novae with the shortest orbital periods T 
Pyx, IM Nor, and Cl Aql have decline time-scales fa significantly 
longer than all the other systems. In sectionwe argue that there 
is some evidence for the other observed properties of these systems 
to be in line with the other recurrent novae. We discuss the mass 
of their ejecta in sections|2d]and|2^ and the duration of their soft 
X-rays emission phases and the velocity of the ejecta in section [23] 
In section]^ we put forward a tentative explanation for this phe¬ 
nomenon. Finally, in section|^we summarize our findings. 


* E-mail: andrea.caleo@astro.ox.ac.uk 


The latest outburst of the recurrent nova T Pyx occurred in 2011 
and the development of its light curve and spectrum has been fol¬ 
lowed more thoroughly than for any other recurrent novae to date 
(e.g. see |Shore et al.|20T2 and Shore et al.|2013| for the optical and 
UV spectra, De Gennaro Aquino et al.||2014 and|Chomiuk et 


2014 


et al. 


for the X-rays, Nelson et al.|2014 for the radio, and Patterson 


2014|for the time-evolution of the orbital period). This allows 


astronomers to compare the features of T Pyx with those of other 
classical and recurrent novae. 


2.1 The mass of the ejecta 

It is difficult to accurately determine the mass of the ejecta of a 
nova. It has usually been estimated based on photoionization mod¬ 
els alone, which requires panchromatic observations (at least in¬ 
cluding the ultraviolet and optical resonance and principal recom¬ 
bination lines). Recently, a new technique has been applied to sev¬ 
eral bright novae including T Pyx, to obtain density, filling factor, 
and mass estimates based on electron density determinations us¬ 
ing forbidden lines from the nebular stage spectra (see [Shore et al.| 
|2012| and |Shore et al.|2013y In brief, either or both of the isoelec- 
tronic forbidden lines of [N II] and [O III] provide determinations 
of the density and electron temperature. The electron density map 
is obtained by line profile ratios instead of using integrated fluxes 
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2 Andrea Caleo and Steven N. Shore 


Table 1. The Galactic recurrent novae from|Schaefer|p010| with the addition of the latest eruptions of U Sco (2010) and T Pyx (2011). The systems are listed 
in order of orbital period. We have not included the eruption of T Pyx close to 1866, inferred by|Schaefer, Pagnotta & SharalpOlO). 


RN 

Lpeak (mag) 

K,i„ (mag) 

*3 (d) 

Sorb (d) 

Eruption yeai's 

T Pyx 

6.4 

15.5 

62 

0.076 

1890, 1902, 1920, 1944, 1967, 2011 

IMNor 

8.5 

18.3 

80 

0.102 

1920, 2002 

Cl Aql 

9.0 

16.7 

32 

0.62 

1917, 1941,2000 

V2487 Oph 

9.5 

17.3 

8 

~ 1 

1900, 1998 

U Sco 

7.5 

17.6 

2.6 

1.23 

1863, 1906, 1917, 1936, 1945, 1969, 1979, 1987, 1999, 2010 

V394 CrA 

7.2 

18.4 

5.2 

1.52 

1949, 1987 

TCrB 

2.5 

9.8 

6 

228 

1866, 1946 

RS Oph 

4.8 

11 

14 

457 

1898, 1907, 1933, 1945, 1958, 1967, 1985, 2006 

V745 Sco 

9.4 

18.6 

9 

510 

1937, 1989 

V3890 Sgr 

8.1 

15.5 

14 

519.7 

1962, 1990 


assuming that the nebular spectra are (roughly) isothermal. The line 
profiles are modeled adopting a bipolar structure with a linear ve¬ 
locity law jRibeiro et al.|2013) l to obtain the geometry of the isother¬ 
mal ejecta. Knowledge of the distance from the object allows us to 
convert the H/3 and Ha fluxes to emission measures that are then 
used to obtain the filling factor using the independently obtained 
electron densities and volumes. This is then used to obtain the total 
mass of the ejecta. The temporal development of the electron den¬ 
sity can also be determined using multi-epoch modeling ( [Schwarz] 
[ 20 ^ . 

Studies have been conducted for Cl Aql by |Iijima]]2012^ and 
for T Pyx by [Shore et al.| ( (20T^ , resulting in the value Mej Ri 2 • 
10~® Mq for both systems. This number is on the lower end of the 
scale for ejecta masses of classical novae. [Shore et al^]2013^ also 
quote a value for the filling factor of T Pyx, / « 3 • 10“^. To our 
knowledge, there is no estimate of the mass of the ejecta of IM Nor. 

We note that the value of the mass of the ejecta of T Pyx is 
controversial. While the spectroscopic value is low, there are indi¬ 
cations that support a higher value, the most suggestive of which is 
based on the orbital period change undergone by the system during 
the eruption jPatterson et al.|2014[ ( and gives a mass of the ejecta 
Mej > 3 • 10“^ Mq. We argue that the dynamics of the binary 
system during the eruption is complex and would deserve further 
consideration, discussing this issue in appendix A. 

Another indication of a possible high value for the mass of 
the ejecta of T Pyx is provided by radio observations ( [Nelson et al.[ 
[2014^ . The peculiarity of T Pyx as a nova is evident in this spectral 
range, as the radio flux began rising surprisingly late (~ 50 days 
after the outburst), and no simple model of instantaneous, homol¬ 
ogous explosion is able to satisfyingly fit the data and provide an 
estimate for Mej. However, a model designed to fit to the later part 
of the light curve gives Mej « 4 • 10“^ Mq, and more elaborate 
models for the whole light curve suggest a range of (1 — 30) • 10“® 
Mq. 


2.2 Mej — t 2 correlation 

Low-mass ejecta have lower density than high-mass ones and they 
are also expanding at higher speed so that their density will drop 
faster. The ionization of the ejecta will therefore be fast if the mass 
of the ejecta is low, and the time for the decline will be accord¬ 
ingly short, so that a correlation between these quantities can be 
expected. Despite the uncertainties in the masses, [Della Valle et al.[ 
\2QQ2) found a correlation using a sample of 18 novae with more 
than one independent source for the mass of the ejecta of most sys¬ 
tems. These results are shown in figure [7|[Della Valle et al.[(|2002[l 


oo 

o 



Figure 1. Optical decline time-scale t 2 and mass of the ejecta for a set of 
novae; M^j is measured in units of 10“® Mg and t 2 in days. The approxi¬ 
mate positions of the recurrent systems T Pyx, Cl Aql, and U Sco, based on 
estimates of t 2 from the light curve of T Pyx, on the light curves reported 
by|SchaeferH2010j and on the masses of the ejecta estimated byjShore et al.| 
]2(il3|,|Iijima|j2012| and [Diaz et al.H2012| (Mej > 3 ■ 10“® Mq for U 
Sco) have been marked in red. V382 Vel is highlighted in the article by[Della[ 
[Valle et ar]j2002[ because that system was one of the main objects of study 
of their article. The dashed curve is given hy equation ^ and it is the best 
linear fit for the Log(Mej) — Log(t 2 ) relation. 


provide the fitting relation (with a 95% confidence level): 

Log(Mej) = 0.274(±0.197) • Logjfa) -b 0.645(±0.283), (1) 

with Mej measured in 10“® Mq and t 2 in days. 

Applying equation Q with Mej = 2-10 ® Mq gives a very 
short predicted decline time-scale: the higher value in the interval 
permitted by equation Q is f 2 < 1 d, while the observed decline 
time-scales can be estimated from the light curves as f 2 ~ 20 d for 
Cl Aql and f 2 « 30 d for T Pyx. This places Cl Aql and T Pyx 
significantly off the curve of equation 0. We conclude that the 
optical decline time-scales of these systems are significantly longer 
than those of systems with a similar value of Mq. 

2.3 X-Rays emission 

Classical novae, including T Pyx, IM Nor, Cl Aql, and other recur¬ 
rent systems, show an extended strong soft X-rays emission after 
the outburst which is modeled as continuing nuclear reactions in 
the white dwarf envelope. The peak photon energy is Speak < 1 
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Figure 2. Correlation between the velocity of the ejecta and the duration of 
the soft X-rays phase. The x-axis shows the FWHM velocity of the Ha and 
iT/3 spectral lines near visual maximum, an indication of the ejecta velocity. 
The y-axis gives the time at which the ratio of soft over hard emission drops 
below a certain threshold and is an indication of the turnoff time of the nu¬ 
clear reactions at the white dwarf surface. Filled circles are known recurrent 
novae. Half filled circles are suspected recun'ent novae based on their char¬ 
acteristics. The arrows indicate upper (pointing down) and lower (pointing 
up) boundaries. The dashed line shows the fit by [Greiner, Orio & Schar- 
|teT|j2003| based on a subset of the repoiled data. Adapted from [Schwarz 
<201 !> . 


keV (hence the often used term supersoft), and the FUV/X-rays 
spec trum resembles a blackbody distribu tio n only at low res olution 
(see Orio, Covington & 6gelmanj2Q01 


and 


Schwarz|201 1 1 . Since 


the ejecta of the system are also expelled from the white dwarf at 
the beginning of the eruption, a correlation is expected between the 
duration of the supersoft emission, the mass of the ejecta, and the 
decline time-scale of the system. 

[Schwarz] ( |2011^ performed a correlation study of the velocity 
of the ejecta, which is related to their mass and can be directly ob¬ 
tained spectroscopically, and the duration of the X-rays emission. 
We show their results in figure]^ where we have marked the posi¬ 
tion of T Pyx (see section [?.3.1| l. Cl Aql, and IM Nor. In the sense of 
this correlation, the supersoft turn-off times of these three systems 
are not especially anomalous. 

|Schwarz| ( |2011^ have also studied the correlation between the 
decline time-scale t 2 and the turn off of the X-rays emission. Their 
results are shown in figure]^ where we have again marked the po¬ 
sition of T Pyx, Cl Aql, and IM Nor. Although the scatter is signifi¬ 
cant and there is a scarcity of systems with long decline time-scales, 
some correlation is evident. T Pyx, IM Nor, and Cl Aql are clearly 
deviant in this correlation. We conclude from the analysis of fig¬ 
ures]^ and that while the X-rays emission and ejecta velocities of 
these systems are similar to those of the other recurrent novae, their 
decline time-scales t 2 are anomalous. 


2.3.1 A word of caution: emission from T Pyx 

The X-rays emission from T Pyx has been studied with the Swift and 
Suzaku facilities. The results are discussed by jPe Gennaro Aquino] 
jet al.| ( |20T4l l and jChomiuk et al.|p014| >. The X-rays behaviour of T 
Pyx is anomalous, including the time evolution of the hard-to-soft 
X-rays emission ratio which is used by jSchwM^ ( |201 1) to define 
the turn off time of the phase of soft X-rays emission. |De Gen-j 
jnaro Aquino et al^ <2014|l suggest that the emission received from 
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Figure 3. Correlation between the optical decline time-scale t 2 and the 
turnoff time. As in figure]^ filled circles are known recurrent novae; half 
filled circles are suspected recurrent novae based on their characteristics; 
the arrows indicate upper (pointing down) and lower (pointing up) bound¬ 
aries. The approximate position of T Pyx is depicted as a red circle, showing 
its peculiar behaviour. The upper boundaries for Cl Aql and IM Nor, already 
included in the analysi s by[Schwai^j2Ql l) , have been highlighted. Figure 
adapted from [Schwarz| l |201 1| . 


T Pyx might have been produced by the hot ejecta and not by the 
central source, while [Chomiuk et al.[ ( |2Q14) discuss the possibility 
of a dual-phase ejection, with a two-months interval between the 
ejections. 

To include T Pyx in the correlation studies of the other no¬ 
vae, we adopted the turnoff time toff ~ 140 d, again given by |De| 
[Gennaro Aquino et aI.| ( |20T4l l. The previous arguments would still 
hold for the other two systems if T Pyx were to be excluded for its 
peculiarity. 


3 REPROCESSING OE THE RADIATION DURING THE 
OUTBURST 

In light of these observations, we propose the following scenario 
for the peculiar behaviour of T Pyx, IM Nor, and Cl Aql: 

• Mass transfer from the companion star onto the white dwarf 
resumes immediately after the eruption with transfer rate M. 

• A fraction / of the transferred material orbits in the Roche 
lobe of the white dwarf, with density distribution p(r, t). 

• In most systems, after a time At that is possibly as short as a 
few days, radiation from the white dwarf ionizes the ejecta and the 
optical luminosity declines. 

• In a few systems, whose particular properties we will discuss, 
the accreted material of mass fMAt absorbs the radiation through 
photoionization before it hits the ejecta, impeding their ionization. 
The energy is then reemitted at a much lower temperature, and is 
less effective in ionizing the ejecta. 

We illustrate this scenario schematically in figure]^ which shows 
the state of the system at the start of the optical emission (on the 
left) and after a time At (on the right). For convenience, we have 
not drawn the figure to scale: the radius of the ejecta a few days 
after the eruption in a scaled figure would be considerably larger. 

This screening effect is favoured in systems with a high mass 
transfer rate, that implies higher density within the Roche lobe of 
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the white dwarf. The mass transfer rates of T Pyx and Cl Aql are 
known to be high. The co nventionally accepted values are Mt Pyx = 
(1-^5)- Mq yr~^ dSelvelli et al.|2008| and MciAqi « 10"'^ 
Mq yr“^ jHachisu, Kato & Schaefer|2003t. There are no estimates 


Godon et al. 


1 2014 1 , based 


of M for IM Nor. A recent analysis by 
on the fit of disc models to the data of far ultraviolet spectroscopy, 
proposes even higher values of Mrpyx. of order 10“® Mq yr~^. 

We propose that this effect is also favoured for systems with a 
low orbital period for the following reason. A low orbital period is 
related to a smal separation, and a small volume of the Roche lobe 
of the white dwarf. Thus, for a given accumulated mass, the region 
around the white dwarf should be at higher density p{r,t). Since 
the effectiveness of the orbiting material in absorbing the radiation 
depends on the recombination rate of its ionized atoms, which is 
proportional to , a high density is a key element for the occur¬ 
rence of this screening effect. 

It is very difficult to devise a comprehensive model of the ir¬ 
radiation of the ejecta taking into account the transfer of material 
into the Roche lobe of the white dwarf, because the hydrodynamics 
of the flow of matter during the peak of the irradiation is very com¬ 
plex. Such model is beyond the scope of our present paper. Instead, 
we present a schematic model for estimating the amount of material 
required for the shielding of the ejecta to occur. 


3.1 The model 

We assume that the mass transfer rate after the eruption is the same 
as in the quiescent state (see appendix B for discussion), and that the 
radiation from the white dwarf is time independent and has a black- 
body distribution. We assume that material around the white dwarf 
has a spherically symmetric distribution p{r) and uniform tempera¬ 
ture r,„. Our aim is to determine the minimum value of the fraction 
of material / that is needed to absorb the radiation incoming from 
the white dwarf after a time Af since the eruption. Any result / > 1 
would imply that the effect we are describing is unrealistic. 

There are few a priori constraints for the value of / from hy¬ 
drodynamics. To our knowledge, there are no models of accretion 
disc formation during the peak of the irradiation of a nova. General 
simulations, as those by [Lanzafame, Belvedere & Molten!] ( |2006| l 
(see their figure 2), still allow a great range of values for /, which 
varies from effectively 0 for high viscosity material to about 10~^ 
for very low viscosity values. 

We adopt the following fixed numerical values for the system: 

• For the white dwarf, a mass of 1 Mq, and a radius of 0.01 
Rq. For the secondary, a mass of 0.2 Mq (these are reasonable but 
debatable values for the T Pyx system; see|Uthas, Knigge & Steeghs| 

( |M0l l). 

• Temperature of the white dwarf Twd = 7 ■ 10® K, correspond¬ 
ing to Eddington luminosity. 

• Time that would be required for the ionization of the ejecta 
without the intervening matter Af = 5 d, which is a typical decline 
time-scale of very fast classical novae (|Strope, Schaefer & Henden| 

|20T0) . 

• Temperature of the material in the Roche Lobe of the white 
dwarf: Tm = 10^ K. This value is consistent with energy balance in 
the Roche lobe; see appendix Cl for discussion. 

• For the opacity, we neglect the neglect the metals in the chemi¬ 
cal composition of the transferred material and adopt for the hydro¬ 
gen and helium abundances: X — 0.75, Y = 0.25. This facilitates 
the calculations adopting a lower limit for the total opacity. See ap¬ 
pendix C2 for a discussion. 


We label as N^, (r) the number of photons per unit frequency at 
a distance r from the white dwarf, and proceed to determine (r) 
as the number of photons emitted from the surface of the white 
dwarf minus those absorbed by the material in the Roche lobe. The 
photons are absorbed by processes of photoionization when they 
strike the atoms; we neglect the scattering processes because they 
are approximately elastic and have a smaller effect on N^{r). 

At a temperature Tm = 10^ K, part of the hydrogen is ionized 
because of the collisions and does not contribute to the absorption. 
The amount of ionized atoms depends on the density of the material. 
Our model takes this effect into account by considering a reduced 
number of H atoms, given by the fraction of atoms that are not colli- 
sionally ionized, as derived by solving the Saha equation. We label 
such fraction as F(p). At Tm = 10'^ K, for the densities relevant 
to our scenario, only a negligible fraction of the He atoms are ion¬ 
ized by collisions. We acknowledge that our description constitutes 
an oversimplified view of the complex interaction between the ra¬ 
diation from the white dwarf, the collisions between the atoms and 
those between the atoms and the electrons, and the recombination 
processes, but we proceed in this way as our aim is to obtain an 
order-of-magnitude estimate of / rather than an exact value. 

We label as B^{Tv/d) the black body spectrum at the white 
dwarf surface temperature, as ai{v) the photoionization cross sec¬ 
tion for photons of frequency v striking atoms of species i, as (T) 
the case B recombination coefficient of the species i, as np, Ua and 
rie the number densities of H, He nuclei (i.e. protons and a parti¬ 
cles) and electrons. In what follows, ni(r) is the number density of 
the species i that contribute to the recombination - photoionization 
balance, so that rip = F{p)np and fia — Ua- We denote as Trl 
the radius of the Roche lobe and compute it by means of Eggleton’s 
formula ( |Eggleton|1983| l. 

Analytical formulae for an, a He and for the cross section 
are given in |Draine'| ( [201 l| l. The recombination coeffi¬ 
cient aHe++ is derived from the hydrogen recombination coefficient 
as described in section 4.2 of |Osterbrock & Eerland| ( |2006t . The He 
photoionization cross section (v) is interpolated from the 

data by |Yan, Sadeghpour & Dalgamo|(T998^ . The photoi onization 
cross section of He"*" is given in section 5.4 of 


Cox 


120001 . 


The number densities of protons (either free or in H atoms) 
and a particles (either free or in He atoms) are given by: 


np(r) = Xp{r) —, 

?72p 


( 2 ) 


na{r) = Y p{r) — , (3) 

Tflcy 

where mp and rria are the masses of the H and He nuclei respec¬ 
tively. 

The equations that give the number of photons radiated from 
the white dwarf and their absorption by the material in the Roche 
lobe are: 


Wi,(7?wd) = 47ri?wD^i/(2wD)^, (4) 

= ^ ai{v)ni{r)^, (5) 

i=H,He,He+ 


The equations that give the balance between the recombination and 
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start of the optical emission 
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a few days later 



Ejecta 


^ Suppressed Radiation 



Figure 4. Schematic illustration of the scenario discussed in section for T Pyx, IM Nor, and Cl Aql. The optical emission starts at the end of the fireball 
phase. The material in the Roche lobe of the white dwarf has been blown away by the eruption, but the mass transfer from the companion through the Li point 
is uninterrupted. After a few days, the expansion of the ejecta has caused their density to drop. However, the material in the Roche lobe of the white dwarf 
suffices to absorb most of the radiation (illustrated by an'ows in the figure) and re-emit it at a much lower temperature, decreasing its effectiveness. The ejecta 
ai‘e not ionized and the optical emission continues. The figure is not to scale: the radius of the ejecta a few days after the eruption in a scaled figure would be 
considerably larger. 


the ionization of H, He, He"*" and He^"*" are: 


3.2 Results 


(^m) — 


nn 


'^He+^eQ^He+(^m) — 


riHe 

47rr^ 


+ 00 


^He++^^^He++(*^ni) — 


^ + oo 

Jo ' 

I 

i: 


N„a^^^^^+{v)dv, 




( 6 ) 


(7) 


^He+ 


iV^O-He+^He++(j^)rf*^, (8) 


where the dependence of hjj+ = njj+ (r), njjj+, njjg++, rie and 
on r has been omitted for readability. 

Finally, the equations of conservation of the total number of H 
and He nuclei and electrons are: 

%+ (r) + nuir) = F{p)np{r), (9) 

%e++(f) + «.He+(f) + nHe(r) = Heir), (10) 

( 1 - F(p) )np + nH+ (r) + %,+ (r) + 2nHe++ (r) = n, (r). 

(11) 

Solving equations (5) - (11) in the variables N„, uh, nH+ ’ 
nHe+> h-He++’ with the boundary condition (4), we determine 
(r). We thus obtain /3, the energy that is free to escape the Roche 
lobe: 

P = hvN,,{RpC)dv^ / ( y hvN,,{R-^p,)dv^ ■ (12) 


To solve equations (5) - (11), a density distribution law p{r) is re¬ 
quired, for a total mass in the Roche lobe fMAt. Since the recom¬ 
bination rate depends on the square of the density, the required value 
of / will depend strongly on the density profile. Unfortunately, the 
distribution of the material depends on the hydrodynamics of the 
accretion process and is therefore difficult to predict in advance. 

We discuss the simplest case of uniform density distribution 
p(r) = p. It is important to note that this is the least favourable 
choice since the effectiveness of recombination grows strongly with 
p and the radiation is therefore more easily absorbed if there are 
regions of high density and regions of low density rather than a uni¬ 
form profile. As an illustrative example, we also present results for 
the case p{r) oc . The profile has physical significance 

because it is the density distribution that would result in the case 
of steady, free-fall, spherical accretion onto the white dwarf, and 
represents an extreme case worth studying. 

As stated in section [TT[ we assumed Mwd = 1 Mq, Mcomp = 
0.2 Twd = 7 • 10® K. We extracted the orbital periods from 
table[^and set Mtpv x = 3 ■ 10~^ M(7) vr~ ^, in the mid-point of the 
range of values from 


Selvelli et al. 


I ^08|_^nd iUciAqi — 10 


M,7 


by |Hachisu, Kato & Schaefer j2Q03| . 


yr“^ as given 

We consider the material in the Roche lobe to be effective in 
absorbing the radiated photons if 90% of the emitted energy is ab¬ 
sorbed, i.e. if P = 0.1. We denote as fc the fraction of transferred 
material that has to orbit in the Roche lobe for this to occur and 
give values of fc for T Pyx and Cl Aql in the uniform density and 
p{r) oc cases in tablej^ 
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uniform density 

p(r) oc r 

TPyx 

fc = 2.5 ■ 10“^ 

/c = l.l-10-^ 

Cl Aql 

M 

1 

o 

q 

II 

fc = 2.6 ■ 10“^ 


Table 2. Fraction fc of transferred material that needs to stay in the Roche 
lobe of the white dwarf of T Pyx and Cl Aql in order to absorb 90% of the 
radiated energy. 


TPyx 

Cl Aql 

f 

/3 

f 

0 

1.0 

10“^ 

0.9 

4.0 

10“^ 

0.7 

2.0 

10“2 

0.4 

5.0 

10“2 

0.4 

2.5 

10“2 

0.1 

6.0 

10“2 

O.I 

3.0 

10“2 

0.0 

7.0 

10“2 

0.0 


Table 3. Fraction / of material in the Roche lobe and fraction of non- 
absorbed energy /3 for T Pyx and Cl Aql, in the case of uniform density 
distribution. 


In all cases, a fraction /c < 0.1 is sufficient to absorb the 
outgoing radiation. For p{r) oc /c is of order 10~^. If the 

density distribution differs greatly from the ones we studied with 
our model, even smaller values of / might suffice. Thus, the trans¬ 
ferred material seems to be sufficient to shield the ejecta. In light of 
these results, we argue that the scenario we described is realistic. 

We show in table|^the f — ft correspondence for some values 
around fc for both T Pyx and Cl Aql, in the case of uniform density 
distribution, ft depends strongly on / so that a moderate variation 
in the amount of absorbing material can have a large effect. 

We argued that a longer orbital period P^b should lead to 
higher value of /c. We show some values of Porb — fc in table 1^ 
assuming a mass transfer rate M = 3 • 10“® Mq yr”'^, as for T 
Pyx, and uniform density. As we expected, /c grows with Poi-b- For 
systems with an orbital period of several days, not even all of the 
transferred material would suffice to shield the ejecta. 

The results for ft depend uniquely on the product /M, not on 
the individual values o f f and M. For the values MTPyx ^ 10“ 
Mq yr“^ proposed by 


Godon et al. 


1 2014 1 , lower values of / are 
required to shield the ejecta, with fc of order 10“® for the uniform 
density model and 10“'^ for the p oc model. 


4 CONCLUSION 

In this paper, we discussed some properties of the three recurrent 
novae with the lowest orbital periods, T Pyx, IM Nor, and Cl Aql, 
noting that they have the longest optical decline time-scales but the 
mass and velocity of their ejecta and the duration of their X-rays 
emission are similar to those of systems with more rapid optical 
declines. 


Porb (d) 

fc 

0.05 

or 

O 

1 

to 

0.1 

3.4 ■ 10“2 

0.2 

7.6 ■ 10-2 

0.3 

1.2 ■ 10“i 

0.5 

1.6 ■ 10-1 

1 

3.3 ■ 10-1 

2 

6.7- 10-1 

3 

1.0 


Table 4. Fraction /c of material corresponding to ft = 0.1 for the mass 
transfer rate of T Pyx, as a function of the orbital period, in the case of 
uniform density distribution. 


We put forth a scenario to explain the reason for this occur¬ 
rence. We propose that some of the material transferred from the 
secondary star during the eruption contributes to the absorption of 
the radiation from the white dwarf and shields the ejecta before their 
ionization occurs. Quantitative evaluation of this scenario indicates 
that for a system with a large mass transfer rate and short orbital 
period, the amount of transferred material is indeed sufficient. A 
more thorough analysis of this problem requires understanding the 
hydrodynamic of the transfer process. Unfortunately the situation 
is so complex, not to mention the presence of the ejecta and the 
irradiation, that the feasibility of any model on this subject is ques¬ 
tionable. 

The detailed physics of the ejecta has not been addressed in 
this study, but it plays a key role in determining the photometric 
properties of the system. While it is clear that the radiation repro¬ 
cessed from the material in the Roche lobe of the white dwarf is 
reemitted at a much lower temperature, and is therefore less effec¬ 
tive in ionizing the ejecta, its effect on the ejecta should be studied 
to gain a deeper understanding of the light curve of the system. 

Another key aspect that remains to be investigated is the con¬ 
sequences of the shielding effect on the spectrum of the system. 
Although the spectrum of T Pyx during various phases of its out¬ 
burst has been observed (see |Shore et al.|2012||Shore et al.|2013t , 
the major obstacle in understanding the behaviour of this particular 
system is the complexity of the spectra rather than the scarcity of 
the data. The analysis of the spectra prompted [Shore et al.| ( |2013| ) 
to suggest that T Pyx is peculiar with respect to most of the classi¬ 
cal novae and that its ‘extended opaque phase could be due to the 
formation of a cool common envelope after the explosion, causing 
a recombination wave to move outward through the ejecta that also 
extinguishes the XR emission. It could then slowly clear as the WD 
settles into a stage of quasi-static nuclear burning and develops 
a supersoft source The scenario we described provides an expla¬ 
nation for how such an envelope could be formed and motivates 
further thought on the spectra of the system. 
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APPENDIX A: ORBITAL PERIOD OE T PYX BEFORE 
AND AFTER THE ERUPTION 

[Patterson et al.|p014| l reported a campaign to track the photomet¬ 
ric wave of T Pyx and follow the evolution of its orbital period, 
gathering data over the 1996-2011 period. They reported a gradual 
increase of the orbital period during the quiescent phase and a jump 
of -1-0.0054(6)% at the eruption. They claim that this implies that 
the mass of the ejecta of T Pyx is at least 3 • 10“® Mq. We argue 
here that the variation of the orbital period doesn’t imply a lower 
bound to the ejected mass. 

[Patterson et al.[ ( [^14[ ) state that ‘during the eruption, mass loss 
should increase Porh, and angular-momentum loss should decrease 
it’. A positive change in the period would imply that the mass loss 
effect wins and the mass loss is at least 3 • 10“® Mq. Although no 
derivation is shown, this result appears to be based on the analy¬ 
sis by [Livio[([1991^, which in turn is based on that by[Shara et al.| 
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( |1986| l. However, a tacit but fundamental assumption is at the core 
of the paper by Shara et ah: the formulas reported, most notably 
their equation (5), are only valid if the orbit of the system before 
and after the eruption is circular. While it is possible to imagine a 
gradual circularization of the orbit after the eruption, so that this as¬ 
sumption might hold for the study of the long-term evolution of the 
system, the eccentricity right after the eruption is likely non-zero. 

In general, the variation of the period of a system in a circular 
or elliptical orbit depends on the change of the masses of its com¬ 
ponents and the mechanical energy of the system. Equation (11) of 
|van den Heuvel| \1992) relates the change in the period P of a bi¬ 
nary system to the change of its total mass M and its energy per 
unit reduced mass t: 


dP dM 3 de 
~P ^ ~M ~ 2~T 


(Al) 


We express this result in terms of the masses mi, m 2 of the compo¬ 
nents of the system and in terms of the total energy P = /re (where 
^ — (1/mi -I- l/m 2 )~^ is the reduced mass): 


~P 


dmi (, ^ 3 m 2 \ , dm 2 f, , 3 mi \ 3 dP 

~WV 2m7/ 2m^/ 2 ”^ 


(A2) 


For mass loss from the primary star, dmi = —M^j, dm 2 = 0: 


5P /"i 3m2\ 3 5E 

1^ ^ ~W 2 rrd) ~ 


(A3) 


where 5P is the period variation and SE is the change of mechan¬ 
ical energy of the system. It is very complex to even estimate SE 
from models of the ejection mechanism (during which nuclear en¬ 
ergy is converted to thermal energy through degenerate H burning 
at the surface of the white dwarf, and then to mechanical energy of 
the ejecta and possibly the binary system); even the kinetic energy 
of the ejecta is currently not very well known (see e.g. the discus¬ 
sion in |Shara et al.|2010) . Without an estimate of SE, equation l |A3[ l 
cannot provide an estimate of Mej. 


APPENDIX B: MASS TRANSFER DURING THE 
ERUPTION 

In the analysis of section|^ the mass transfer rate from the compan¬ 
ion to the white dwarf during the first few days of the eruption was 
considered to be the same as that at quiescence. We discuss here 
this assumption. 

The secondary star is strongly irradiated during the outburst. It 
could be argued that this would cause an expansion of the star and 
an increase in the mass transfer rate. The problem of irradiation- 
induced mass transfer in novae has been widely discussed in the 
literature; however, the focus has always been on the long term evo¬ 
lution of the system rather than the first few days after the outburst. 
The most influential work in the field was conducted by |Kovetz, Pri-| 
|alnik & Shara| ( [T98^ . These authors determined the depth reached 
by the irradiation and argue that the main scattering process in the 
atmosphere of the star is electron scattering. They obtained for the 
mass of the irradiated layer ~ 5 • 10“®M, where M is the 
mass of the companion. The companion would then inflate on a 
time-scale of « 0.1 yr, overfilling its Roche lobe and causing a 
significant increase in the mass transfer rate. 

In contrast, we are considering the evolution of the system on 
a very short time-scale, the first few days after the outburst. Since 
the star would not have had time to fully expand, the scenario pro¬ 
posed by |Kovetz, Prialnik & Shara| ( |1988^ does not apply to our 
case. Moreover, we argue that the incoming photons are absorbed 


and re-emitted by processes of photoionization rather than electron 
scattering and give an independent estimate of Min. 

We achieve an order-of-magnitude estimate of M[„ by assum¬ 
ing that the atmosphere of the companion is composed of hydrogen, 
with a density profile similar to that of a non-Roche-lobe-filling 
main sequence star of the same mass (we acknowledge that this is 
inaccurate in a neighbourhood of the Li point). We make use of the 
density profile of a star with mass M ~ 0.2 Mq, effective tem¬ 
perature Teff = 3200 K, gravity Log{g) = 5, and solar chemical 
composition extracted from the NextGen models (see |Hauschildt,| 
[Allard & Baron|(T999| >). We evaluate the effect of a monochromatic 
irradiation at the peak wavelength of the emis sion from the white 
dwarf, at Two = 7 ■ 10® K as in section 3.1 The resulting peak 
photon energy is Speak ~ 60 eV. 

Proceeding as in section [uT| we label as Nj{z) the number 
of photons per unit area and unit time that reach depth 2 ; in the 
atmosphere, as nuiz) and nn+{z) the number densities of H and 
that contribute to the recombination - photoionization balance, 
as (Tph the photoionization cross section of H atoms at Speak ~ 60 
eV, as S(p) the fraction of H atoms that would not be collisionally 
ionized at density p and temperature Tm, and as a{Tm) the type B 
recombination coefficient of hydrogen. As in the case of sectionjUT] 
we assume Tm = 10"^S; see appendix Cl for discussion. We label 
as a the distance between the white dwarf and the secondary star. 
The equations that determine N,y{z) are: 


NAO) 

Z/WD 1 

4:7T(1^ _£/peak 

(Bl) 

dN-y 

dz 

= -N-/{z)nH{z)aph, 

(B2) 

n^i+iz 

)Q;(7|n) ■— ■ '^H (■2)o'ph-/V'Y ( 2 ^) , 

(B3) 

nuiz) 

-|-nH+( 2 ) = np{z), 

(B4) 

np{z) 

TPi' ( \\/^NextGen (-2) 

= i^(pNextGeii(2:)) -• 

mp 

(B5) 


We solve equations (B2) - (B5) with the boundary condition (Bl), 
and define the penetration depth Zi„ as the depth at which 90% of 
the incoming energy has been absorbed. The result is z^ ~ 10^ 
cm. 

Assuming uniform irradiation of the exposed of the secondary 
star, we estimate Mi„ as: 


M„ = 2TrRl [ ” p{z)dz ^ 1O"^'‘M0. (B6) 

Jo 

This result is six orders of magnitude lower than that of| Kovetz, Pri-| 
jalnik & Shara| ( [T988) l. The effect of the irradiation on the secondary 
star is therefore very small and the mass transfer rate is unlikely to 
be affected by it. 


APPENDIX C: ASSUMPTIONS IN THE MODEU OF 
SECTION 3.1 

Cl Temperature in the Roche lohe of the white dwarf 

In section [UT| we argued that the temperature Tm of the material in 
the Roche lobe of the white dwarf is Tm ~ 10^ K. We discuss here 
this assumption. 

We use an energy argument to determine an approximate value 
for Tm- The material in the Roche lobe is heated because of the 
radiation incoming from the white dwarf. It dissipates this energy 
by processes of photo-recombination to excited states that release 
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photons that are less energetic than the ones from the white dwarf. 
The mean free path of these photons in the Roche lobe is: 


I 


mfp — 


1 

pt^bf(P; Tin) 


(Cl) 


We use Kramer’s formula for the bound-free opacity Kbf, assuming 
metallicity Z = 0.02 and H abundance X = 0.75. For the case of 
uniform density in the Roche lobe of T Pyx described in section [T^ 
with f ~ fc = 2.5 • 10“^, this gives = 3 • 10® cm. The radius 
of the Roche lobe is of order trl ~ lO'^® cm. This means that, 
depending on the position in the Roche lobe, a significant fraction 
of the re-emitted photons are not re-absorbed and contribute to the 
cooling of the system; most of them, however, are not lost, so that 
our scenario is only approximately correct. However, Kbf(p, Tm) oc 
Tn^®'® depends strongly on Tm, so that, for a temperature T = 1.2 ■ 
10"^ K, (mfp and r*RL are of the same order. To obtain a temperature 
estimate, we assume transparency for photons from recombinations 
and that these are the main source of cooling. 

The energy per unit time and volume absorbed from the in¬ 
coming radiation is: 


1 

£abs(r) = -—/ N^hvdu. (C2) 

47r7’ dv Jq 

The energy per unit time and volume emitted by the material 
is: 


£em = A(r)nenH+ , 


(C3) 


where A{T) is the cooling function. We use the values tabulated by 
Sutherland & Dopita ^1993f for solar composition (i.e. [Fe/H] = 


0). Such cooling function was computed under conditions valid for 
the interstellar medium rather than a stellar environment with con¬ 
siderably higher density (of order prl ~ 10“'^® g cm~® in the 
Roche lobe of the white dwarf); however, in both cases the main 
contribution to the cooling function at T > 10^ K is from hydro¬ 
gen recombination to excited states (see e.g. |Dalgamo & McCray| 
(TTTII i), so that the predicted cooling rate is still approximately 
correct. Figure 8 of [Sutherland & Dopita| ( |T993| l shows the cool¬ 
ing function for a large range of temperatures and various chemical 
compositions. 

We solve the equation £cm{r,T) = £abs(i') in the variable T 
for every r with the values of ne(r) for the uni¬ 

form density case for T Pyx with / = /c. This temperature profile 
is not self-consistent, as Eabs has been obtained by assuming a uni¬ 
form Tm = K profile, and will not reproduce in details the 
real profile in the Roche lobe. The temperature maximum in figure 
IClI is located where the bulk of the radiation is absorbed from the 
material in the Roche lobe; here a significant number of atoms are 
neutral and the cooling is most efficient. 

Our results show that the temperature in the Roche lobe, 
though not uniform, does not deviate significantly from 10^ K. This 
is due to the fact that the cooling function A(r) varies very steeply 


with T for T > 10"^ K (see again figure 8 of Sutherland & Dopita 


1993 i: a temperature profile that spans a range of a few 10^ K is 


suited to describe regions with significantly different rates of emit¬ 
ted energy. On the other hand, the recombination coefficients Oi 
and the fraction of non collisionally ionized H atoms F{p) do not 
depend on T in a similar manner, and they can be safely computed 
at the uniform temperature T,„ = 10^ K. This is why we have not 
considered the effects of a temperature gradient in the Roche lobe. 
Solving equations (5) - (11) of section ED with the higher value 
T = 1.2 • 10"^ K gives very similar results for the value of /c and 
for the temperature profile of figure [CT] 



Figure Cl. Temperature profile obtained imposing £cm(r, T) = eabs(r) for 
the case of uniform density for T Pyx with f = fc- 


Although we have focused the discussion of this appendix on 
just the case of section]^ a similar argument holds for the assump¬ 
tion Tm = 10^ K of appendix B. In that case, almost all of the 
photons originated in recombination to excited states contribute to 
the cooling of the system, because the depth Zin ~ 10^ cm reached 
by the irradiation is at low optical depth, as can be seen from the 
data in the NextGen model (r < 10“"^). We conclude that for both 
problems, Tm ~ 10^ K is an acceptable approximation and the tem¬ 
perature gradient is not very significant. 


C2 Metals in the Roche lohe of the white dwarf 

In section ED we neglected the absorption of photons originated 
at the surface of the white dwarf by elements heavier than He. We 
discuss here this assumption. 

Inspection of the solutions of equations (5) - (11) in the /? = 
0.1 cases described in section \T2\ shows that, in all of these cases, 
most of the atoms, whether of H or He, are ionized. In this situation, 
the number of photons per unit volume and time absorbed for effect 
of the species X, equal to the number of recombinations that take 
place, is approximated by: 

dN-y I 

=nx+ncax, (C4) 

at \x 

where ax is the recombination coefficient of the species X. The 
relevance of the species X to the absorption depends therefore on 
the value of the product nx+ ax ■ 

Assuming a metallicity of order 10~®, the total number density 
of the heavy elements with atomic number Z > 6 is of order 10“® 
times that of the H atoms. Some recombination coefficients ax 
have been tabulated, including those of the most common hydro- 
genic atoms ( [Storey & Hummer|1995^ . They are generally higher 
than the H recombination coefficient qh, but most often well within 
a factor of 10® from it. We conclude that the contribution of the 
heavy elements to the absorption is of order 10“^, or lower, com¬ 
pared to that of the H atoms. It is therefore possible to neglect it in 
a schematic model as the one of section |3T| 
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